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ABSTRACT 

We examine the observational consequences of partial gaps being opened by planets in protoplanetary disks. 
We model the disk using a static a-disk model with detailed radiative transfer, parametrizing the shape and size 
of the partially cleared gaps based on the results of hydrodynamic simulations. Shadowing and illumination 
by stellar irradiation at the surface of the gap leads to increased contrast as the gap trough is deepened by 
shadowing and cooling and the far gap wall is puffed up by illumination and heating. In calculating observables, 
we find that multiple scattering is important and derive an approximation to include these effects. A gap 
produced by a 200 (70 M^) planet at 10 AU can lower/raise the midplane temperature of the disk by up to 
~ -25 / +29% ('^ -1 1/ + 19) by shadowing in the gap trough and illumination on the far shoulder of the gap. At 
the distance of Taurus, this gap would be resolvable with ~ 0.01'' angular resolution. The gap contrast is most 
significant in scattered light and at thermal continuum wavelengths characteristic of the surface temperature, 
reducing or raising the surface brightness by up to order of magnitude. Since gaps sizes are correlated to planet 
mass, this is a promising way of finding and determining the masses of planets embedded in protoplanetary 
disks. 

Subject headings: planet-disk interactions - protoplanetary disks — planets and satellites: detection — radia- 
tive transfer 



1. INTRODUCTION 

With the advent of adaptive optics systems on large ground- 
based optical and near-IR telescopes, we are beginning to 
be a ble to image pro toplanetary disks around young stars 
(e.g. iThalmann et al.l 12010. LkCa 15). These gas-rich disks 
are where we expect giant planets to form, since gas giants 
need that large reservoir of gas from which to accrete their 
massive envelopes. However, the interpretation of structure 
in these imaged disks can be problematic, because scattered 
light traces only the optically thin and diffuse surface lay- 
ers of the disks rather than the overall structure of the disk 
(I Jang-Condell & BossI 120071) . On the other hand, if growing 
planets do indeed significantly perturb the disks in which they 
are born, observation of these signatures can mean the detec- 
ti on of nascent planets in protoplanetary disks. 

lJang-Condelll(l2008L henceforth JC08) calculated the shape 
of a dimple created by the gravitational potential of a planet 
embedded in a disk and the accompanying thermal pertur- 
bations in shadowed and illuminated regions of the dimple. 
pang-Condell (2009, henceforth JC09) predicted the observ- 
ability of these dimples. Neither of the papers considered the 
clearing of an annular gap in the disk by tidal forces. Plan- 
ets above about 30 are able to open partial gaps in disks, 
as demonstrated in numerical hydrodynamic simulations (e.g. 
iBate etal.ll2003h . These gaps are much larger in scale than 
the local dimples modeled in JC08 and JC09 and are there- 
fore a promising way of detecting and characterizing planets 
in protoplanetary disks. If the gap is well resolved enough 
to determine its depth and width, then we can determine the 
mass of the planet to within a factors of a few. 

JC09 considered only single scattering and direct thermal 
emission to calculate observable signatures of disk pertur- 
bations. However, multiple scattering, particularly at high 



albedos, can significantly increase disk brightnesses. While 
Monte Carlo methods directly calculate multiple scattering, 
they are computationally intensive so that iteratively cal- 
culating the disk structure is prohibitively time consum- 
ing (e.g. Walker et al. 2004; Dullemond & Dominik 2004; 
Tannirkulam et al. 2008; Pi nte et al.1 120081: Mulders et al.] 
2010). The approach here is to exploit a few analytic approx- 
imations to the solutions of the radiative transfer equations to 
efficiently calculate the disk structure. 

Numerous people have carried out hydrodynamic 
simulations of gap openi n g in disks by planets (e.g . 
Paardekoop er & Papaloizoul |2QQ8|; lAvliffe & Bat^ 120091; 
de Val-Borro et al. I I2006t I Bate et al. [20031: lEdgar & OuillenI 
2008) However, the effects of illuminatio n of these gaps by 
the central star has not been well-studied. I Wolf et alT (120021) 
examine this problem in a flat disk with constant H/r = 0.05 
and find that the far edge of the gap is not illuminated. On 
the other hand, Var niere et al.l (I2006D find that the far edge 
in a flared disk, which is more typical of T Tauri disks, is 
heated and puffed up, creating a positive feedback loop that 
enhances the appearance of the gap. In both these cases, the 
gap was created by a Jupiter-mass planet, which clears nearly 
all the material from its orbital, creating a deep and wide gap. 
In this work, we focus on partially cleared gaps created by 
sub- Jupiter mass planets. 

A full three-dimensional hydrodynamic simulation includ- 
ing radiative transfer is very computationally intensive, so we 
simplify the problem by assuming a fixed disk structure in 
terms of its radial surface density profile, 5](r). To further 
simplify the computation, we assume an axisymmetric disk 
structure. This means that we do not include the behavior of 
spiral density waves, but only the large scale gap opened by 
the planet. 
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2. THE JC MODEL 

The methods used for calculating radiative transfer in pro- 
toplanetary disks for both radiative heating and observable 
quantities are those used in JC08 and JC09. In this work, 
we now apply these methods to an axisymmetric disk struc- 
ture rather than an azimuthally limited local perturbation, and 
we account for multiple scattering in the observables. We will 
refer to these methods as the JC model, to constrast to Monte 
Carlo methods. 

The disk model incorporates radiative transfer based 
on the methods of iJang-Condell & Sasselovl (|2003h and 
iJang-Condell & Sasselovl (I2QQ4 and adopted to self- 
consistently determine the disks thermal and pressure 
structure as described in JC08. Stellar irradiation is an 
important heating source in protoplanetary disks, and the 
amount of heating is sensitive to the angle of incidence at 
the disk surface. If the shape of the surface of the disk 
changes on scales smaller than the disk thickness, the plane- 
parallel approximation fails. To account for this kind of disk 
structure, radiative heating is integrated piecewise over the 
surface of the disk rather than assuming a thin plane-parallel 
disk model. This method of calculating radiative transfer 
is fully-three dimensional and not directionally dependent 
on a fixed coordinate system. More detailed explanation of 
the method can be found in JC08 and references therein. 
Calculation of the observables is done as detailed in JC09, 
with modifications to approximate the effects of multiple 
scattering, as discussed below. In this section, we summarize 
the essential features of the model. 

2.1. Initial Conditions 

We follow the same method as used in JC08 and JC09 to 
calculate the initial disk structure. The stellar parameters are 
mass = IM©, radius = 2.6Rq, and effective temper- 
ature Teff = 4280 K, consistent with a protostar with an age 
of 1 Myr (Sie ss et al.. .2QQQ). We assume a constant-a disk 
model where the vis cosity is parameterized sls u = a^^CsH 
(TShakura & Sunyaevl[T973|) . The disk parameters are accre- 
tion rate M = 10~^Mq yr~^ and viscosity parameter o^ss = 0.01, 
parameters typical for T Tauri stars. 

The initial conditions are calculated in a two-step process. 
In both steps, we iteratively calculate the density and temper- 
ature structure of the disk, including radiative transfer of the 
stellar irradiation at the disk surface. Viscous heating is in- 
cluded, but is only important inwards of a few AU. Beyond 
this distance, the primary heating source is stellar irradiation, 
which depends sensitively on the angle of incidence at the 
surface of the disk. The surface is defined to be where the 
Planck- averaged optical depth to stellar irradiation is 2/3. 

In step 1 , we generate a locally plane-parallel model for the 
entire disk. For this initial locally plane-parallel model, we 
calculate the disk at logarithmically spaced intervals of V2 
from 0.25 AU to 256 AU. 

In step 2, we select a radially and azimuthally limited 
wedge of this disk to calculate in more detail, including the 
three dimensional curvature of the wedge, but assuming ax- 
isymmetry. The model presented here is a slice with radial 
range (r) from 3 to 20.8 AU and vertical distance from the 
midplane (z) from to 5.5 AU in order to model the gap cre- 
ated by at planet at 10 AU. We sample the initial density and 
temperature at 100 grid points in both r and z directions, us- 
ing cylindrical coordinates. In order to ensure that that this 
initial disk is in thermal and hydrostatic equilibrium, we re- 



calculate the heating from stellar irradiation, now removing 
the assumption that the disk is locally plane parallel and in- 
cluding the full three-dimensional curvature of the disk. The 
resulting density profile is So(^)- The total disk mass within 
the simulation boundaries is 0.0025 Mq . 

2.2. Gap Modeling 

The depth of a gap that can be opened by a planet depends 
not only on the mass ratio between planet and star, but also 
on the disk properties. In a low viscosity disk, the critical 
threshold for gap-opening is generally considered to be the 
mass at which the planet's Hill radius equals the thermal scale 
height of the disk. The Hill radius is rnm = (q/3y^^a where 
q = Mp/M^ is the mass ratio of the planet to the star. The 

thermal scale height is H = Cs/^k where Cs = ^/kf~/~i~L is the 
thermal sound speed and = yGM^Jc? is the Keplerian or- 
bital angular speed. Here, k is Boltzmann's constant; T is the 
disk's midplane temperature; /i is the mean molecular weight 
of of the gas, assumed to be primarily molecular hydrogen; 
and G is Newton's gravitational constant. 

The gap opening criterion for a viscous disk has been found 
empirically to be 



where Reynolds number is defined as IZ = r^Q.p/v 
(ICrida et al.l2006l). The ir results are also consistent with those 
of Var niere et aP (12004 . 

For the purposes of this paper, we adopt gap profiles sim- 
ilar to those calculated by lBate et al.l (l2003l) . who model gap 
opening by planets of varying masses in a disk using three- 
dimensional hydrodynamic simulations, focusing specifically 
on planets that only partially open gaps in disks. A gap is 
modeled as an ad hoc perturbation imposed on the initial con- 
ditions. The surface density of a disk modified by a gap of 
width w, depth J, and position a is 

S(r) = Eo(r) {1 - Jexp[-(r-^)V(2w2)]} . (2) 

This is consistent with the results of (ICrida & Morbideilil 
|2007|) . who find that planets whose masses are less than that 
of the disk open local gaps in disks rather than clearing in- 
ner cavities. Since the study presented here addresses planets 
under a Jupiter mass (M/), they are well under the disk mass. 

The depth and width of the ga p compared to pla net mass are 
determined from the results of Bate et"an (|2003|) . In Figure 
[H we reproduce the surface density profiles of gap opening 
by planets from Figure 2 of Bate et al. (2003), having used 
DataThiefl to acquire the values. The unperturbed density 
profile is plotted as a solid line, and goes as So(r) oc r~^/^. 
The planet is located at r/a = 1. The gaps opened by planets 
of 0.03, 0.1, 0.3 and 1 Mj are plotted as dashed, dotted, dot- 
dashed, and long-dashed black lines, respectively. We fit all 
but the 1 M/ planet to Equation ([J]), and plot these fits as ma- 
genta lines in Figure [TJ The best fit parameters are tabulated 
in Table [T] The maximum deviation between the empirical 
gap profile and the fit, expressed as a precentage of the unper- 
turbed disk density, is 0.04% for the 0.03 Mj planet, 3% for 
the 0.1 Mj planet, and 9% for the 0.3 Mj planet. The devi- 
ation for the 0.3 Mj planet is greatest dXr/a < 1, indicating 
asymmetry of the gap that is not accounted for in the model 

^ B. Tummers , |http ://datathief . org | 
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assuming that /Mq = 3x10 




Figure 1. Profiles of gaps carved by planets in disks and the Gaussian fits 
to them. The unperturbed disk profile is plotted as a solid line, while the 
dashed, dotted, dot-dashed, and long-dashed lines are gaps opened by planets 
of g = 3 X l£;^J_x lO^f^x 10""^ , and 1 X 10"^ , respectively. The black lines 
are from IBate et all ( 120031) . and magenta lines are Gaussian fits, following 
Eq. (|2}. The gap opened by the g = 1 x 10~^ planet is not well fit by a a 
Gaussian. 

presented here. Nevertheless, fitting the gap to a Gaussian is 
a useful model for parametrizing the disk response to a planet 
without the need for running a full hydrodynamic simulation. 

In dBate et al.1 12QQ3I) . 7^ = 10^ and H/a = 0.05. Using 
Eq. ([!]) as the gap opening criterion, this gives a gap-opening 
threshold of ^crit = 1.06 x 10~^, or slightly more than 1 Mj. 
For comparison, rnm = H when q = 3.75x 10~^, so the viscous 
gap opening criterion gives a mass more than twice as large as 
for an inviscid disk. In the disk model adopted in this paper, 
at 10 AU the midplane temperature is 50K, H/r = 0.048, and 
Q/gg = 0.01, so 7^ = 4.4 X 10"^ and the gap-opening threshold is 
^^rit = 1.97 X 10"^ or almost 2 Mj. 

Given the difference in disk properties, we cannot assume 
that the same mass planets open the same size gaps in each 
disk. We estimate the masses of planet according to the fol- 
lowing procedure. Defining 



3 H 50 

4 rmii q/Z 



(3) 



we assume that G is the relevant scale for determining the 
depth of the gap opened by a planet, so that planets with the 
same value of G open similarly sized gaps. Then , the 0.03, 
0.1, an d 0.3 M/ planets in the disk modeled by Bate et al.l 
(I2003h have mass ratios of ^ = 3 x 10"^ 1 x 10"^ and 3 x 10"^, 
respectively, and have G values as tabulated in Table [T] For 
the disk parameters adopted in this paper, the equivalent 
planet masses are then approximately 20, 70 and 200 M^, 



Table 1 

Best-fit parameters for gaps opened by planets. 





d 


w/a 


max. 
error^ 


G 


derived q 


derived 
planet mass 


3 X 10-^ 


0.014 


0.078 


0.04% 


18 


6.7 X 10-5 


22 Me 


1 X 10-4 


0.56 


0.11 


3% 


6.2 


2.2 X 10-4 


72 Me 


3 X 10-"^ 


0.84 


0.17 


9% 


2.5 


6.2 X 10-4 


210 Me 


1 X 10-3 


4 


4 


4 


1.0 


1.9 X 10-3 


620 Me 



^As simulated in lBate et al.l ( |2003|) 
bError= (S-Sfit)/So 

^Actual masses used for this work. In the text, the masses have been 
rounded to 20, 70, and 200 Me for convenience. 

^The gap opened by the q= 10~^ planet in the IB ate et al.l ( 12003 ) simulation 
is not well-modeled by a Gaussian. 




Figure 2. Cartoon diagram of radiative feedback on disk structure. The star 
is represented as a yellow disk, and the planet by a blue dot. The disk surface 
represents a contour of constant density. The left image shows the initial gap 
opened in the disk, with the dotted line showing the original, unperturbed 
disk surface. Stellar illumination on the surface of the gap creates shadowed 
and brightened regions. Shadowing and cooling occurs in the disk trough, 
and the far side of the gap is illuminated and heated. The right image shows 
the response of the gaseous disk material to the cooling and heating: the 
shadowed region contracts and deepens the gap, while the illuminated far 
side expands and is elevated. 

Holding I](r) constant and assuming axisymmetry, we re- 
calculate the temperature structure of the disk, using the same 
iterative radiative transfer method as described in ^2.11 The 
effect of the gap is to create a shadow within the trough in 
density created by the gap. The far side of the gap, which 
is now exposed to more direct stellar illumination, is bright- 
ened. A schematic of this is shown in Figure [2l This results 
in cooling within the trough and heating on the far wall. The 
gas that composes the bulk of the disk material responds to 
this cooling and heating by contracting and expanding. This 
changes the vertical density profile of the gap region, and the 
illuminated surface now must be recalculated. For this rea- 
son, the heating and cooling of the disk must be calculated 
self-consistently with the density structure of the disk. 

2.3. Observables 

Calculation of synthetic images of disks with gaps is done 
as in JC09, with the addition of including multiple scattering. 
In our comparisons to Monte Carlo radiative transfer codes, 
we find that multiple scattering turns out to be an important 
component of disk emission. Multiple scattering generally 
has the property of increasing the brightness of a diffuse disk. 
In the case of high albedos, the additional brightening can be 
significant. Here, we present an estimate for the brightness 
due to secondary scattering of both scattered stellar irradia- 
tion and thermal emission from the disk. We assume isotropic 
scattering in these derivations. 

2.3.1. Multiple Scattering of Stellar Irradiation 

Here we consider pure scattering of stellar photons. From 
JC09, the intensity of singly scattered light of frequency v 
from the disk surface is 



(4) 



where ujjy is the albedo, Bj^(T^) is the intensity emitted at the 
stellar surface (here assumed to be a blackbody), r is the dis- 
tance from the star, /i is the cosine of the angle of incidence 
of stellar light, and r] is the angle between the line of sight to 
the observer and normal to the surface. The angle r] is dis- 
tinct from the inclination angle of the overall disk, which we 
represent as /. The angles are illustrated in Figure O 

For multiply scattered photons, we first need to determine 
the diffuse scattered radiation field within the disk, Jg. This 
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Figure 3. Illustration of how inclination affects brightness of scattered light. 
The path length through the illuminated layer represents the relative bright- 
ness seen by each observer. The cosine of the angle of incidence of stellar 
light at the surface is /i. The angle between the surface normal and the ob- 
server is r\. 



was solved for in I Jang-Condell & Sasselovl (I2QQ3L 12004 us- 
ing plane-parallel equations for radiative transfer as 



4/2 l-3gV 
2 + 3/i 



1+2^/3 



exp(-gr^) - 3/i exp(-r^ / ji) 



(5) 



where g = V3(l-a;^) (see alsol D' Alessio et alJ (Il998h ). 

Given J^, then the brightness of photons scattered to the 
observer after two or more scatterings is 



Jo 



QXp(-Tohs)Jsdrohs 



(6) 



where Tobs is the line-of- sight optical depth. If r± is the optical 
depth perpendicularly below the surface of the disk, and r] is 
the angle of the observer with respect to the disk normal, then 
'^obs = ^±/ COS 7^. Then, 



P- 



4/2 1- 
2 + 3/i 



3/i 



(1 +2^/3) (1 +^C0S7y) (l+C0S7^//i) 



(7) 



and 



TS_JS , ,s_ ^l^f^^lBl^(T*) 

-il +i2 - " ^ 



4r2(/i-Fcos7/) 

(2 + 3/i)(/i + cos7/) 



(l + 2g/3)(l+gcos7?) 



}■ 



(8) 



2.3.2. Scattered Thermal Emission 



Assuming a plane-parallel disk atmosphere, the radiative 
transfer equation is 



cos 7/7— =Ij, 



(9) 



where r] is the viewing angle with respect to the line perpen- 
dicular to the disk, and Ty is the total extinction perpendicular 
to the surface, including both absorption and scattering. For 
scattering and absorption, the source function is 



Sy=UJyJy + {\-UJy)By 



(10) 



where ujy is the wavelength-dependent albedo and By = By(T) 
is the local thermal emission. Adopting the Eddington ap- 
proximation. 



3 Sr2 



= (l-iJy)(Jy-By). 



(11) 



We have By from the temperature structure of the disk, and 
we can integrate Eq.[TT]to get Jy. Finally, we integrate Eq.[9] 
to get the emitted intensity at the surface of the disk. 

In the case of an isothermal slab of finite optical thickness 
(Tmax), we Can set the boundary conditions at r = and r = 

'^max 



1 dJ 

J^(r = 0)=-= — (r = 0) 
V3 9r 



(12) 



1 dJ 

V3dr 



■Tmax) (13) 
(iRvbicki & Lightman|[T979h . Solving, 

Jy=By(T){l (14) 

exp [(r-rmax)V3(l -cj^)] +exp [-rV3(l -ujy)] \ 
(1 - V^-ojy) exp [-rmaxV3(l-^i.)] + 1 + V(l-ojjy) J * 

If we assume that T varies slowly with r, then we can ap- 
proxim ate Jy by Eq. ([T4b locally, as done in ID' Alessio et al.l 
(1200 Ih Then the emergent intensity can be found by integrat- 
ing Eq. Q, assuming that the background intensity is 0. To 
simplify the calculation, we assume that the disk is locally 
plane parallel, so r] = /, the inclination angle. For a face on 
disk, r] = 0. 

The total intensity is then ly =1^+11, although as a general 
rule /2 dominates in optical to mid-IR, and II dominates at 
longer wavelengths. 

3. VALIDATION 

The JC radiative transfer model for calculating the disk 
structure for this work is fully three-dimensional, in the sense 
that radiation impinging on the surface of the disk is allowed 
to propagate in all directions throughout the disk. How- 
ever, it relies on a locally one-dimensional analytic solu- 
tion to the radiative transfer equation and adopts a number 
of simplifying assumptions such as m ean opacities (see also 
iJang-Condell & Sasselovl I2003L 120041) . The novel approach 
presented in this paper is the iterative calculation of the self- 
consistent density and temperature structure. As shown by 
the significant temperature perturbations created by shadow- 
ing and illumination on gaps, the self-consistency is an im- 
portant aspect to consider in the analysis of radiative transfer 
in disks. 

In order to validate the radiative transfer prescription 
adopted in this paper, we compare our results to a Monte 
Carlo radi ative transfer calcu lation on the final disk density 
structure (iTurner et al.l 1201 ll). The M onte Carlo approach 
is similar to that of iPinte etaH (120061) . We find the radia- 
tive equilibrium temperatures by following a large number 
of photon packets from the star through scattering, absorp- 
tion and re-emission until escape to infinity. In this way the 
energy is conserved exactly. We su m the r adiation energy ab- 
sorbed all along the packet paths (iLucyl [T999i) and relax to 
equilibrium by choosing each re-emitted packet's frequency 
to adjust the local radiation field for the updated temperature 
(.Bjorkman & Wood .2001.) . The scattering is assumed to be 
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isotropic. Each of the calculations shown involves 10^ photon 
packets. 

In the following section, as we discuss our results for the 
disk temperature structure and observable quantities, we com- 
pare the JC model to the results found by the Monte Carlo 
(MC) method. 

4. RESULTS 
4.1. Gap Structure 




Figure 4. Surface density (top) and midplane temperature (bottom) profiles 
for a disk with and without a gap created by a planet at 10 AU. Solid, dotted, 
and dashed lines indicate planet masses of 0, 70, and 200 , respectively. 
The thermal perturbation is caused by shadowing and illumination by stellar 
irradiation at the disk surface. The black lines indicate the results of the JC 
models while the gray lines show JC models excluding viscous heating at 
the midplane. The magenta lines on the temperature plot show the results of 
a Monte Carlo radiative transfer calculation on the same three-dimensional 
density structure. 

The surface density profiles imposed on the disk for 0, 70, 
and 200 planets are shown in the upper plot of Figure (H 
as solid, dotted, and dashed lines, respectively. The calculated 
midplane temperatures for the corresponding disk models are 
summarized in the lower plot of Figure HI in black for the JC 
models, and in magenta for the MC calculations. 

The Monte Carlo calculations include only the disk annulus 
between cylindrical radii 3 and 20.8 AU. For computational 
speed, the highly optically-thick material within 3 AU is omit- 
ted. Without the obscuring material, the wall facing the star 
at 3 AU would become unexpectedly hot. We mitigate this 
by immediately discarding any Monte Carlo photon packet 
absorbed or scattered for the first time in the innermost col- 
umn of grid cells. Some residual excess heating comes from 
the fact that omitting the inner disk reduces the column to 
the star even in the atmosphere, where the innermost cells are 
optically-thin. This residual excess is significant only inside 
about 4 AU. On the other hand, truncating the disk model at 
the outer radius reduces temperatures, since photon packets 
reaching 20.8 AU can escape to infinity, radiatively cooling 
the midplane. Considering these approximations together, we 
believe the Monte Carlo results between 4 and 18 AU accu- 
rately represent the solution that would be obtained at greater 
computational expense in a disk model continuing inward to 
the sublimation radius and outward to the disk's far edge. 



The JC temperature is consistently higher than the MC tem- 
peratures overall. One cause of this discrepancy is that the JC 
model includes viscous heating at the midplane, so that the 
final temperature is 



t = (t;+t:) 



4x1/4 



(15) 



where Tr and Ty are the temperatures resulting from stellar ir- 
radiation and viscous heating, respectively. The viscous tem- 
perature is given by 



9 GM,M 



327r (jgr^ 



1-1^ 

r 



1/2" 



(r^ + 2/3) (16) 



where is the Stefan-Boltzmann constant and is the ver- 
tically integrated optical depth, rd = -J^ XRp^z, with respect 
to Xr^ th^ Rosseland mean opacity. 

In Figure IH we plot Tr in grey to show the effect of re- 
moving viscous heating. Viscous heating contributes more at 
smaller r, and is neglible beyond ~ 15 AU. Interior to 5 AU, Tr 
is less than that calculated by the MC model, and outside of 5 
AU, the MC temperatures fall more rapidly than the JC tem- 
peratures. Thus, viscous heating can explain the discrepant 
temperatures only at small radii. Another possible explana- 
tion is that the difference lies in the treatment of opacities in 
the two models, as will be described below. 

The gap opened by the 20 Mq planet is less than 2% in 
depth, and produces no more than 2% excursions in temper- 
ature from the unperturbed disk, so small as to be negligi- 
ble. We thus set 20 Mq as the lower bound on a detectable 
planet in this disk at 10 AU, or more generally, planets with 
q < 0.036 ^crit or G > 18 are undetectable, whereas planets 
with q > 0.12^crit or G < 6 do significantly perturb the disk. 

In contrast to the results of JC08, the midplane tempera- 
ture is significantly affected by gap opening. JC08 considered 
planets up to 50 in the absence of a gap, but here we find 
that a 70 planet should open a significant gap. In a disk 
with lower viscosity, a 50 might open a similarly sized 
gap. This shows the importance of including large scale non- 
linear dynamical interactions in planet-disk models. 

Figure m illustrates that a gap in a disk creates an S-shaped 
perturbation to the midplane temperature profile. The tem- 
perature at the position of the planet is lowered, although the 
temperature minimum is inward of the planet position. This 
is because cooling at the gap minimum is mitigated by heat- 
ing on the far gap wall. In the JC model, the minimum and 
maximum temperature deviations at the midplane for the 70 
Me planet are -6 K (-11%) at 8.9 AU and +8 K (+19%) at 
13.4 AU, respectively. In the MC model, these values are -6 
K (-13%) at 8.0 AU and +9 K (+25%) at 13.1 AU, respec- 
tively. For the 200 planet, the values are -13 K (-25%) 
at 9.1 AU and +12 K (+29%) at 13.9 AU in the JC model, 
and -11 K (-23%) at 7.9 AU and +14 K (+38%) at 13.4 AU. 
In general, the magnitude of heating and cooling in Kelvins 
is comparable, but the percentages differ because the unper- 
turbed temperature of the MC model is smaller. The positions 
of the minima and maxima in the JC model are slightly inward 
of the minima and maxima in the MC model, which might be 
caused by the same reasons that give rise to the lower mid- 
plane temperatures in the MC model. 

The full radial and vertical temperature structure of our 
models are shown in Figure \5\ The vertical axis is stretched 
compared to the horizontal axis in these plots. The top, mid- 
dle, and bottom plots show gaps created by 0, 70, and 200 
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M0 planets, respectively. The JC models are plotted in black, 
while the MC models are plotted in magenta. 





10 15 
radius(AU) 



Figure 5. Temperature cross-sections of disk models with and without gaps. 
Black and magenta lines show the JC and MC methods, respectively, of cal- 
culating radiative heating. Top: A disk without a planet or a gap. Middle: 
A disk with a gap created by a 70 planet. Bottom: A disk with a gap 
created by a 200 planet. 

The abrupt shifts in temperature in the upper left corners 
of the JC model are purely a numerical artifact. This region 
is extremely low density, so this region is considered to be 
outside the bounds of the simulation, and the temperature is 
assigned based on zero optical depth and a radial dependence. 
Since these regions are very low density, the temperature dif- 
ferences have no effect on hydrostatic equilibrium nor on the 
observables, and do not affect the remainder of the analysis 
presented in this paper, though the topic should be addressed 
in future studies. 

The disk surface is coincident with regions where the tem- 
perature changes rapidly in the vertical direction, i.e. where 
the temperature contours lie close together. Above the sur- 
face, the JC model is generally warmer than the MC model. 
However, at or just below the surface, the MC model is gen- 
erally warmer. Closer to the midplane, the JC model again 
becomes warmer, as discussed previously. The warmer tem- 
peratures in the MC model below the disk surface indicate 
that heating by stellar irradiaton penetrates deeper in the MC 



model than the JC model. The probable cause of this and 
other differences in temperatures between models, is that the 
JC model calculates the propagation of stellar photons into the 
disk using mean opacities, which treats all stellar photons as if 
they have the same opacity. In the MC models, opacity varies 
with wavelength, so photons with wavelengths longward of 
the stellar peak have a lower opacity and can heat the disk 
to deeper depths. The inclusion of viscous heating in the JC 
model leads to the increase in temperature at the midplane. 

When a gap is imposed on the disk, as in the lower two 
plots of Figure [51 shadowing and cooling within the gap re- 
sulting in the temperature contour lines dipping downward. 
The depth of these dips are comparable in both the JC and MC 
models, although the exact shape differs. This difference may 
again be explained by the the choice of mean opacities ver- 
sus wavelength-dependent opacities. Allowing longer wave- 
length photons to penetrate further into the disk would have 
the effect of smearing out the temperature differences casued 
by surface effects. 

A close look at the surface temperature contours at r > 
20AU reveals that the surface temperatures are lower for the a 
gapped disk in both the JC and MC models. This is a result of 
shadowing of the outermost radii by the puffed up outer edge 
of the gap. This shadowing and cooling effect is quite real, but 
since the shadow extends beyond the simulation boundaries, 
quantifying this effect further is beyond the scope of this pa- 
per. A similar shadowin g by the il luminated outer edge of a 
gap created by Jupiter by Turner e t al. ( 201 IJ. 

4.2. Simulated Images 
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Figure 6. Radial surface brightness profiles of face-on disks with and with- 
out gaps at the wavelengths imaged in Figures [7] & [H as indicated. The 
gapless disk surface brightness are plotted with solid lines. Surface bright- 
ness profiles of gaps created by 70 and 200 planet are plotted as dotted 
and dashed lines, respectively. Black lines indicate the JC models includ- 
ing multiple scattering, while gray lines show only single scattering or di- 
rect emission. For comparison, results from MC calculations on the identical 
disk density structure are overplotted as magenta lines. The 1 mm brightness 
profiles is noisy because of low photon statistics. MC calculations give con- 
sistently higher brightnesses because of multiple scattering and temperature 
differences, as discussed in the text. 
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The radial surface brightness profiles for each of the disks 
are shown in Figure [6l at wavelengths of 1 - 1000 /im. The 
JC models are plotted in black, and the Monte Carlo models 
are plotted in magenta. The brightness profiles of the unper- 
turbed disk are plotted as solid lines, and gaps created by 70 
and 200 planets are plotted as dotted and dashed lines, 
respectively. We show general good agreement between the 
models. The noise in the brightness profiles at 1000 microns 
in the MC models is caused by low photon statistics. This is 
because the 1000 microns wavelength lies out in the tail of the 
thermal distribution. Most photon packets have wavelengths 
near the peak emission of star or disk. 

For comparison, we plot as gray lines in Figure [6] the JC 
brightness profiles for single scattering and direct thermal 
emission only, as calculated in JC09. Inclusion of multi- 
ple scattering of both scattered stellar photons and thermally 
emitted photons bring the JC brightness profiles in line with 
the MC ones. At 1 micron, the albedo is ujy = 0.91 and 
the disk emission is effectively purely scattered light, so the 
disk about twice as bright. At 10 microns, = 0.45, and 
scattered photons and thermally emitted photons are roughly 
equal. Because of the relatively low albedo, multiple scat- 
tering increases the brightness by about 30% averaged over 
radius. At 30 microns, ujjy = 0.45 also, but because thermal 
emission is dominant at the wavelenght, the overall brightness 
is increased by about 70%. Multiple scattering is more effec- 
tive at increasing thermal emission because it allows photons 
emitted in optically thick regions to percolate up through the 
disk. Scattered light, on the other hand, is limited to the pho- 
tons incident on the surface of the disk. As the albedo in- 
crease, the amount of brightening increases: at 100, 300, and 
1000 microns, the cj^ = 0.57,0.82, and 0.95, and the bright- 
ness is increased by factors of 2.6, 4.9, and 9.8, respectively, 
as averaged over radius. 

The perturbation from the gap shadow is evident across all 
wavelengths, reflecting the cooling and heating that take place 
in the shadowed trough and brightened rim of the gap. The 
puffing up of the far rim of the gap is significant enough to 
shadow and cool the outer disk, as evident in the steeping of 
the brightness profile toward 20 AU at 1-100 microns. The 
brightening of the far rim of the gap i n the 1 //m b rightness 
profile is consistent with the results of IVarniere et al. (2006), 
with feedback from cooling and shadowing resulting in a puff- 
ing up of the gap edge. 

Despite the offset in the midplane temperature maxima and 
minima, the appearance of the gap at all wavelengths is con- 
sistent. The JC model tends to overpredict the depth of the 
gap, particularly at 1 micron, and from 30-300 microns. At 1 
micron, this maybe because photons can scatter from the back 
wall of the gap into the gap, and this is not captured in the 
JC model. From 30-300 microns, the difference may arise 
from the detailed differences in the gap temperature structure 
between the JC and MC models. At least qualitatively the 
models agree. 

Having validated the JC model by comparison to the MC 
model, the remainder of our analysis will focus on the JC 
model. 

In Figures |7] and [H we show synthetic images of face-on 
disks without and with gaps at wavelengths from 1 /im to 1 
mm, as indicated. The star is omitted in each image, but its 
brightness is that of of a 4280 K blackbody with radius 2.6 Rq. 
Although the apparent brightness varies as the inverse square 
of the distance, the surface brightness in Jy/asec^ is distance- 



independent. 

The top row in Figure [7] shows the scattered light images 
of the disks at 1 /im. The 3 /im image is very similar to this 
one, just scaled to the stellar brightness at 3 /im and modified 
by the albedo according to Eq. ([8]). The stellar brightness at 
1 /im at a distance of 140 pc is 0.79 Jy. The shadow in the 
gap and the brightening outward of the gap are quite appar- 
ent. The puffing-up of the outer edge of the gap also creates 
a further shadowing of the disk beyond the gap. This disk 
self- shadowing is also evident at longer wavelengths, result- 
ing from cooling in the shadow. Because the disk is not well- 
modeled outside 20 AU, it is unclear whether or not the full 
extent of the disk beyond the gap rim is shadowed or not. 

The middle row in Figure [7] shows simulated images of the 
disks at 10 /im, where the stellar brightness is 0.055 Jy at a 
distance of 140 pc. In the inner, warm region of the disk, the 
brightness is dominated by thermal emission from the disk 
surface. The surface brightness profile is steepest at 10 /im 
because the outer disk is too cold to emit efficiently. At about 
10 AU and beyond, the brightness profile becomes dominated 
by scattered light. Both shadowing of scattered light and cool- 
ing of the disk surface contribute to the imaged gap structure. 

The bottom row in Figure [7] shows the thermal emission 
from the surface of the disk at 30 /im. The stellar brightness at 
this wavelength at a distance of 140 pc is 0.0068 Jy. As with 
the 10 /im image, the radial falloff in surface brightness re- 
flects the radial temperature gradient at the disk surface. The 
gap contrast is highest at 30 microns. This is because the sur- 
face temperatures of the disk are ~ 100 K, where the black- 
body peak is 30 microns. Thus, observations at 30 microns 
will be the most sensitive to temperature perturbations at the 
surface of the disk. At shorter wavelengths, the gap contrast 
is driven by shadowing and illumination at the surface rather 
than actual temperature changes. 

In Figure m the top, middle and bottom rows show thermal 
continuum emission at 0.1, 0.3, and 1 mm, respectively. At 
these wavelengths, the stellar emission is negligible. The disk 
becomes more optically thin toward longer wavelengths, so 
the different wavelengths probe different layers in the disk 
as well as different temperature regimes. In contrast to the 
dimple images calculated in JC09, the gaps are still apparent 
at 1 mm wavelengths because the shadowing effect in the gap 
is on a much larger scale than the localized dimple. Thus, 
searching for gaps in disks with ALMA is a promising way of 
finding planets in during the planet formation epoch. 

One might expect that the shadow produced by the gap is 
offset by the brightened outer edge of the gap, creating a net 
zero effect on the spectral energy distribution (SED). Because 
the simulated region is limited in radius, a full SED cannot 
be produced for the entire disk, as any emission from interior 
to 4 AU or exterior to 20 AU will be omitted. Nevertheless, 
we can assess how a gap in a disk changes the SED by in- 
tegrating over the face on disk images. We approximate the 
contribution from the inner 4 AU by assuming that the bright- 
ness profile goes as , is normalized to the brightness at 4 
AU for the gapless disk, and is unaffected by the presence of 
the gap. The star's emission is approximated as a blackbody 
of 4280 K and radius 2.6 Rq. 

In Figure [51 we show the resulting SEDs assuming that the 
system is at a distance of 140 pc. Since the radiative transfer 
calculation includes only the annulus lying 3 to 20 AU from 
the star, the flux approaches that of the whole disk only at 
wavelengths from 10 to 100 microns, where the annulus emits 
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Figure 7. Simulated images of disks from JC models with and without gaps at 1 (top), 10 (middle), and 30 (bottom) microns. The left, center, and right images 
in each row show disks with 0, 70, and 200 planets, respectively. The contours are spaced at intervals of factors of two in brightness, and the colors trace 
brightness according to the displayed color bars. 
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Figure 8. Same as Figure|7] but at 0.1 (top), 0.3 (middle), and 1 (bottom) mm. 
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Figure 9. SEDs of face-on disks with and without gaps from JC models. 
The SEDs are calculated by summing the total thermal emission over the 
simulated disks. Top: the disk integrated spectrum of a disk without a gap 
is plotted as a solid line, while the disks with gaps carved by 70 and 200 
M0 planets are plotted as dotted and dashed lines, respectively. The dot- 
dashed line shows the stellar photosphere. Bottom: the difference between 
the gapped disk models and the gap-less disk are plotted in magnitudes. The 
70 and 200 models are shown as dotted and dashed lines, respectively. 
The maximum deviation from the gapless model occurs at around ~ 40 /im, 
corresponding to the peak of emission for temperature of the disk surface at 
the gap radius. 

most strongly. Nevertheless, the difference in fluxes should be 
representative. The upper panel shows that the gapped SEDs 
(dotted/dashed line for 70/200 M^) are generally brighter than 
the gapless SED (solid line), particularly at around 30-40 /im, 
corresponding to the peak in thermal emission from the disk 
surface at the gap radius. 

The lower panel of Figure [9] shows the difference between 
gapped and gapless SEDs in magnitudes, with negative values 
being brighter. These magnitude differences should be inter- 
preted as upper limits on the amount of brightening that can 
occur, because emission from the outer disk beyond 20 AU 
has been omitted, and this disk contribution may swamp out 
the relative differences. Shortward of 10 /im, stellar flux (dot- 
dashed line) dominates the SED so the disk contribution is 
negligible. Beyond the submillimeter regime, emission from 
the outer disk becomes increasingly important, so the effect 
of the gap is relatively less than that shown. At best, a 70 
(200) M0 planet increases the brightness of the disk by 0.15 
(0.3) mag at 44 (40) jam, or by 15 (32) %. Beyond about 0.7 
mm, the Amag becomes positive for the 200 gap, indi- 
cating that there is more dimming than brightening at radio 
wavelengths. Since the maximum grain size assumed for cal- 
culation of the opacities was 1 mm, this is likely because at 
wavelengths longer than the maximum grain size in disks, the 
brightness is more representative of the dust distribution in the 
disk rather than the temperature. 

In general, the conclusion is that the brightened outer rim of 
the gap more than compensates for the lack of emission in the 
gap trough. The reason for this is the greater amount of heat- 
ing that occurs on the outer shoulder of the gap, as described 
in M.ll Moreover, the contribution to the brightening is am- 
plified because the surface area of the outer rim is greater than 
the surface area of the gap shadow. However, the contribution 
of the gap to the SED is subtle, and not easily distinguished 



from a more massive disk. The SED alone cannot be used to 
identify a partially cleared gap. 

5. DISCUSSION 

In JC08, we carried out similar calculations on local per- 
turbations to the disk caused by less massive planets than the 
gap producing ones modeled in the present work. Since gaps 
are much larger in physical scale, they produce both a more 
significant temperature perturbation on the disk and a larger 
footprint on the disk itself. Thus, planets that produce the 
gaps modeled in this paper are very likely to be observable, 
whereas the much more subtle perturbations modeled in JC08 
are much harder to observe. Imaging is necessary for detect- 
ing these gaps, since the effects are too subtle to identify in 
the SED alone, as described in ^4.21 

The width of the gaps modeled here are on the order of a 
few AU. If we assume that we need a resolution equivialent 
to 1-2 AU, then if the disk is at the distance of Taurus, 140 pc, 
then the angular resolution required would be 0'^01, indepen- 
dent of wavelength. Scattered light observations face the chal- 
lenge of suppressing the starlight to a small enough working 
angle to observe the gap, in addition to the angular resolution 
challenge. ALMA should have sufficient angular resolution to 
resolve these gaps, but the question them becomes whether or 
not it will have sufficient sensitivity given the long baselines 
and sparse u-v coverage required. The size of the gap scales 
with distance, so a gap at 100 AU would be resolvable with 
O'M resolution. 

We have only modeled gaps at 10 AU in this present work, 
but we can make some extrapolations to how gaps at other dis- 
tances in the disk would behave. The same planet mass opens 
a smaller gap at larger radii because the disk scale height in- 
creases faster than the Hill radius of the planet. Thus, the tem- 
perature variations would be relatively smaller as the planet is 
moved further out. On the other hand, gaps created by plan- 
ets at larger radii should be more easily detectable because 
gap widths scale with distance, relaxing the angular resolu- 
tion requirement, and the inner working angle is larger. For 
interferometers, sensitivity also improves with smaller base- 
lines. We found that the greatest gap contrast was seen at the 
blackbody peak associated with the temperature of the sur- 
face of the disk. At 10 AU, the surface temperature is around 
100 K, so the wavelength of greatest contrast is 30 microns. 
The surface temperature is driven by stellar irradiation, so it 
should vary roughly as r~^/^. This means that to some extent, 
the wavelength of observation may be tuned to the planet's 
orbital radius. 

6. CONCLUSIONS 

We have calculated the thermal effects of gap-opening by 
planets in protoplanetary disks. The thermal feedback leads to 
depression of gap troughs and puffing of gap walls, enhancing 
the observability of gaps carved by planets forming in disks. 
Planets less than one tenth of the critical gap-opening mass, 
or 10 M0 at 10 AU, do not create significant gaps in disks. 
However, a modest gap of only 50%, created by a planet ~ 
10% of the viscous gap-opening mass or 70 at 10 AU, can 
induce a significant perturbation to the temperature profile of 
a protoplanetary disk. These gaps are observable in scattered 
light and thermal emission. 

The ability to determine the masses of planets in disks, to- 
gether with the age of the disks and the locations of gaps, puts 
vital observational constraints on the time scales of planet for- 
mation. If we find that massive planets form early, this might 
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indicate that giant planets form via gravitational instability, 
which is much faster than the competing paradigm of core ac- 
cretion. On the other hand, if we only see gap-forming planets 
at late ages, this might indicate that core accretion is the dom- 
inant process. Upper limits on the luminosity of a planet em- 
bedded in the disk will determine whether giant planet form 
more like brown dwarfs with a hot start (Baraffe et al. 20Q3h 
or more quiescently with a cold start ( Marley et al. 2007V 

Temperature variations in the disk produced by shadow- 
ing and illumination can have profound effects on the form- 
ing planets. We address the consequences for Type I mi- 
gration in a future paper. These temperature perturba- 
tions can also affect the condensation and sublimation of 
volatiles. Planetesimals just interior to the gap may become 
enriched in frozen volatiles, while those just outside the gap 
might become depleted in volatiles via the cold finger effect 
(IStevenson & Luninell988|) . This can affect the overall distri- 
bution of volatiles in the protoplanetary disk, shift the snow 
line, and lead to enhanced planet formation in the volatile en- 
riched regions. 

Topics for future work include carrying out radiative trans- 
fer modeling in a similar way on three-dimensional hydro- 
dynamic simulations of gap-clearing by planets rather than 
relying on a simple analytic model. This would allow us to 
capture non-axisymmetric gap characteristics, particularly the 
region immediately around the planet itself. Our methods can 
also be applied to the inner walls of fully-cleared inner cav- 
ities in transitional disks. If the inner walls are sufficiently 
puffed up, they can shadow the outer disk and create flat radial 
surfa ce brightness profiles seen in some disks ( Grady et al.l 
120051) . Polarized intensity imag es are a promising way t o 
resolve protoplanetary disks (e.g. lOppenheimer et al.l 12008*). 
but geometrical effects make these images difficult to inter- 
pret (e.g. Jang^CondelT^Kuc^^ 

In a forthcoming paper, we will address how inclination and 
polarization affects images of disks with and without gaps. 
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